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Abstract 

We present an overview of the Density Matrix Renormalization Group and its con- 
nections to Quantum Groups, Matrix Products and Conformal Field Theory. We 
emphasize some common formal structures in all these theories. We also propose 
two-dimensional extensions of the variational matrix product ansatzs. 

I. INTRODUCTION 

The Density Matrix Renormalization Group (DMRG) is a powerful numerical Real Space 
RG method introduced in 1992 by S.R. White which can be applied to a large variety of 
Quantum Lattice Hamiltonians defined in Id, quasi Id (ladders) and large clusters |]]. The 
DMRG was originally proposed in the domain of Condensed Matter Physics, where it has 
already become a standard method specially for Id systems, but its range of applicability 
has been extended also to Statistical Mechanics, polymers, Chemical Physics, etc... (for a 
review see ]2|]) and one may expect that it will be applied in the near future to other domains 
of Physics as Quantum Field Theory, Nuclear Physics, etc. 

There are by now excelent reviews on the DMRG || and other related real space RG 
methods @], J5J so it is not the purpose of the present contribution to duplicate material 
already present in the literature. Instead we shall try to give an overview of the DMRG 
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method in order to explore its relationships with Group Theory and its quantum defor- 
mation, Conformal Field Theory (CFT) and the Matrix Product Method. The relation 
between the DMRG and Quantum Groups was suggested in || by studying the RG method 
of q— group invariant Hamiltonians. We also review the variational approach of the DMRG 
in terms of the so called matrix product (MP) ansatzs introduced by Ostlund and Rommer 
@ and suggest a 2D versions of it which may lead to a 2D formulation of the DMRG. 

REAL SPACE RG METHODS: GENERALITIES 

Let us suppose we have a discrete system with N sites and that at each site there are 
two possible states, say spin up and down for a spin system, or occupied and unoccupied 
for a spinless fermion. The dimension of the Hilbert space will grow as 2 N , which makes 
very hard the study of the large N limit, unless some special trick is used. The RG method 
provides a general systematic approach to handle problems with a large number of degrees 
of freedom on the basic assumption that only a small number of states is needed in order 
to describe the long distance physics. How to choose the most representative degrees of 
freedom out of a miriad of states is the central issue of the RG method. This can be done 
in several manners. We shall introduce below some of them and establish their comparison. 

The DMRG method was originally formulated as a real space RG method although 
it admits also a momentum space formulation ||. We next introduce the basic concepts 
common to any real space RG method and later on we confine ourselves to the DMRG. 

The real space RG consist essentially of 3 steps: i) blocking, ii) truncation and iii) 
iteration. First one divides the system into blocks, then one finds an effective description of 
these blocks in terms of intra-block and inter-block interactions and finally one iterates the 
algorithm. One can distingish between 3 types of blockings Kadanof-blocking, Wilsonian- 
blocking and DMRG-blocking. 
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Kadanof Blocking 



We shall first consider the case of a linear chain with N sites. In the first step of the 
Kadanof blocking one divides the system into blocks of 2 sites, thus for a N = 16 chain one 

gets, 

(• •) (• •) (• •) (• •) (• •) (• •) (• •) (• •) (1) 

If every site describes two states then the block (• •) describes 4 states. Eq.(^) is nothing 
but a change from a 1-site basis to a 2-sites basis and hence • • is entirely equivalent to 
(• •). The goal of eq.(l) is to prepare the road to perform the first RG truncation. Indeed 
out of the 4 states we may already want to keep a smaller number, say 3, 2 or even 1. We 
can represent symbolically this operation as follows, 

(• •) — > (• •)' (truncation) (2) 

From a formal point of view the blocking operation is captured by putting parenthesis 
(. . .) around the sites subjected to the blocking while the truncation operation is represented 
by ' acting on the corresponding block. Combining eqs. ([I]) and (0) the chain with N sites 
become after the first blocking and truncation, 

(• •)' (• •)' (• •)' (• •)' (• •)' (• •)' (• •)' (• •)' (3) 

If only one state is kept in (f|) then the RG process ends up since there is already a 
single state to represent the ground state of the system. A dimerized spin chain is a typical 
example of this type of states, where (• •)' is given by the singlet formed by two spin 1/2 
of the chain. 

In general however one keeps more than one state per block (• •)' and so one can continue 
the RG method choosing (• •)' the new site •', i.e. 

(• •)' - (4) 
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The renormalized chain has therefore N/2 effective sites •', which can again be blocked 
as in (|TJ), 

((• •)' (• •)')' ((• •)' (• •)')' ((• •)' (• •)')' ((• •)' (• •)')' (5) 

Performing two more blockings and truncation operations we finally get 

((((• •)' (• •)')' ((• •)' (• •)')')' (((• •)' (• •)')' ((• •)' (• •)')')')' (6) 

Eq. (|^) is the final step of the RG method since the whole chain with N sites has been 
reduced to a single effective site whose dynamics can a priori be easily found by solving the 
final renormalized effective Hamiltonian. 

Wilsonian Blocking 

In his solution of the Kondo impurity problem Wilson || introduced a numerical RG 
method where a single block is grown by adding momentum shells following the so called 
onion scheme. The real space version of this method is summarized in the following eq. 

(((((((((((((((• •)' •)' •)' •)' •)' •)' •)' •)' •)' •)' •)' •)' •)' •)' •)' (7) 

where we have used the same notations as for the Kadanof blocking. While the Kadanof 
blocking follows the pattern B^Bi — > B' 2e , the Wilsonian scheme follows the pattern Bi • — > 
B' g+l i where B g denotes a block with I sites. It seems that the two schemes are completely 
unrelated. Notice however that the number of left and right parenthesis in eqs. (|6]) and 
(0) is the same, namely N — \ = 15. What is different is the order of the brackets. The 
condition for the Kadanof and Wilsonian blockings to be equivalent can be formulated as 
follows, 

((B 1 B 2 )'B 3 ) , = (B 1 (B i B 3 )y (8) 

where Bi(i = 1,2,3) denote generic blocks containing one or more sites. In particular for 
N = 4 one can prove using @ 
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((••)' (••)')' =(((••)' •)'•)' (9) 

Eq. (H) is reminiscent of the associativity of the tensor product of representations in group 
theory and more precisely in quantum group theory ( see below). This equation certainly 
holds if there is no truncation of degrees of freedom, i.e. (B\ B 2 )' = (Bi B 2 ) in which case 
it amounts to the equivalence between different basis. In group theory the relation between 
different basis is given by 6 — j symbols. 

Quantum groups are g-deformations of classical groups ( q = 1 in this notation) where 
some of the commutator and addition rules are deformed (for a review see ||10|| ). The 
representation theory of quantum groups, when the deformation parameter q is generic 
is analogue to that of classical groups. However when q is a root of unit things change 
completely. First of all there are a finite number of regular irreps and the tensor product of 
them is also truncated while keeping the associativity condition (§). 

The existence of an associative truncated tensor is a common feature of the DMRG, 
quantum groups and Conformal Field Theory (CFT) (more on this point below). 

DMRG Blocking 

There are two DMRG algorithms to study open chains. The infinite system algorithm 
uses the superblock Bg • •Bf to grow the chain from both sides according to the Wilsonian 
scheme [p]]. The block Bp* is then truncated to a new block B' e+1 . In this manner the size 
of the system grows indefinitely until one reaches a fixed point beyond which the numerical 
results reproduce themselves. This method is very good in computing bulk properties of the 
system like the ground state energy density. In many cases however it is more convenient 
to study finite size systems whose large distance properties can later on be obtained trough 
finite size scaling techniques. This is notably the case of gapless systems. 

The DMRG algorithm used in these cases is called the finite system method and it is 
extremely accurate. The first steps of this method uses the infinite system algorithm to 
grow both sides of the chain independently until the left and right blocks are a half the size 
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of the chain. The chain with N (even) sites is then obtained by joining a left block B^/ 2 and 
a right block -B^ 2 as follows (the superindex R in B^, 2 indicates that it can be obtained 
from the reflection of a righ block), 

((((((((• •)' •)' •)' •)' •)' •)' •)'(• (• (• (• (• (• (• •)')')')')')')')' (10) 

The superblock (|T0|) is used in the DMRG to enlarge the left block from B N / 2 to B N / 2 +i 
while the right block is reduced from B^/2 to B^, 2 _ u in which case we get 

(((((((((• •)' •)' •)' •)' •)' •)' •)'•)' (• (• (• (• (• (• •) / )') / )') / ) / )' (11) 



If the associativity eq.(||) would hold then the blockings (|10D and (|TTj) would give the 
same result for the GS of the whole chain but of course this is not the case. The next step 
is to again enlarge the left block at the expenses of the right one until one reaches the right 
hand side. There one reverses the trend and grows the right blocks at the expenses of the 
left ones. After several sweeps of this back-and-forth algorithm the GS energy and GS wave 
function converge to a fixed values which are independent of the size of left and right blocks. 

In this moment the splitting of the chain into left and right blocks is independent of their 
size so that the associativity constraint ( ^ is effectively fullfilled. 

The analysis performed so far is rather formal but helps to abstract the blocking proce- 
dure which is common to all the real space RG methods. As a by-product we have shown 
that the blocking and the iteration procedures have to be considered as combined strategies 
to achieve the same goal which is to reduce the whole system to a single effective site. From 
a formal point of view blocking is like tensoring representations. In this sense the RG steps 
can be seen as "putting parenthesis" around the blocks. An exact RG method would be the 
one for which the final result would be independent on the way the parenthesis are put on. 
These lead us to the associativity constraint (|8|), whose fullfillement is the actual goal of any 
exact RG method. 
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THE STANDARD RG ALGORITHM 



In all the real space RG methods there is an algorithm to truncate the collection of two 
blocks BiB 2 down to a new effective block B[ 2 where B\ or B 2 may stand also for a single 
site. 

The standard RG algorithm consist of the following steps 1) diagonalization of the Hamil- 
tonian H BlB2 for the combined block BiB 2 , 2) truncation to the lowest energy states of H BlB2 
and 3) change of basis to the new states kept and renormalization of the old Hamiltonian. 



This method leads to a lot of problems whose origin was first pointed out in reference [11 



following a suggestion by Wilson. Studying the very simple problem of a particle in a box the 
authors of reference Jll| interpreted the bad performance of the standard RG as been due 
to an incorrect treatment of the boundary conditions applied on a block by its neighbours. 
In other words the truncation BiB 2 — > (BiB 2 )' has to take into account the presence of say 
a third block B 3 in contact with the former ones. The key idea is to consider a superblock 
B 1 B 2 B 3 where the effect of B s into the other two blocks can be properly considered. An 
alternative RG-solution to the particle-in-a-box problem, which also takes into account the 
effect of boundary conditions has been given in ||12||. 



THE DMRG ALGORITHM 

Let us choose a superblock made out of three blocks B\B 2 B 3 . The middle block is taken 
to be a single site or two sites B 2 = • or • •. Then one constructs the Hamiltonian H BlB2Ba 
describing the dynamics of the superblock and finds out a given state called the target state, 
which is usually the ground state of the superblock which can be written as 

= J2 V'iiwa |*1,*2,*3> (12) 

where i\, . . . run from 1 up to mi, . . .. The superblock can be regarded either as ((BiB 2 )B 3 ) 
or as (Bi(B 2 B 3 )). Correspondingly the target wave function can be written in two different 
manners, 

7 



J.. _ y. IT. . £(12)3 y . Ho\ 

where U and V are matrices which "diagonalize" the wavefunction and satisfy the orthogo- 
nality conditions, 

Eili2 ^hi2,a ^hn,ot' ~ $a,a' (14) 
Sis V a,i 3 V *',i3 = 5 *,*' (15) 

We have used in (|H|) the singular value decomposition (SVD) of a matrix [|TJ. Z)( 12 ) 3 and 
£)i(23) are s j n g U i ar values of tp regarded as a (m!m 2 ) x m 3 or as a m^mj x m 3 ) matrix. 
Eqs.(|T3|) are the clue of the DMRG method. Let us imagine for a moment that -D^ 12 -* 3 and 
Dp are zero for certain values of a and /3. In this case it is clear that we can truncate the 
states of -B1-B2 ( resp. -B2-B3) down to a smaller set of states a (resp. (3) for which D^ 3 ( 

1 ( 23) 

resp. Dp ) is non zero without loosing any information in order to reconstruct the target 
state ip. Rather than performing the SVD of ip it is more convenient to define the density 
matrices for the subsystems (12) and (23) inside the whole system (123) ill], 

Phil,i' 2 i' 3 = ^ii ^ixiiis ^1*2*3 

Now using (|T3|), (p~5|) we get, 



12) ,,=y [/*. (d^ 3 ) U>> (17) 



(23) _ T/* f n^ 23 )^ T / 

Pi 2i3 ,i' 2 ^ - h,p V /3,i 2 i 3 \ u p ) V M'A 

Eqs.flTTD means that = (D^ 3 ) 2 are the eigenvalues of the density matrix p( 12 ) 

while U is the unitary matrix which diagonalizes p^ 12 -* (similar properties hold for the density 
matrix p( 23 \) Let us call m the number of states kept per block in a DMRG computation. 
This number typically varies between 10 and 1000 depending on the computer resources. 

The DMRG algorithm consists in choosing the m most probable states a, i.e. the states 
with higher value of w a ( let us sort them as (wi > w 2 > w 3 > . . . > w mim , 2 ) This guarantees 
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the best posible representation of the target state if) for every given value of m. Moreover 
the sum P(m) = J2a=i w a °f the probabilities of the m states kept give a reasonable measure 
of the truncation error (recall that tr// 12 ) = J2a w a — 1 an d hence P{m) < 1). 

In many of the Id models studied with the DMRG it turns out that the probability w a 
is concentrated in a few states and that it decays exponentially fast. This implies that with 
small values of m one can achieve a great accuracy in representing the target state. This is 
certainly the case for systems with a finite correlation length |7| Jl3|Jl4j| . For systems with an 
infinite correlation length one has to study finite systems and adjust the number of states 
kept m to the correlation length due the finite size [jy| . 



THE DMRG VERSUS QUANTUM GROUP THEORY AND CFT 

There are certain formal analogies between the DMRG and the theory of quantum groups 
and Conformal Field Theories (CFT) which we shall review below. First of all the DMRG 
truncation of states in (BiB 2 )' has strong similarities with the truncated tensor product of 



irreps of a g-Group where q is a root of unit |16| . Let us choose for example the quantum 



group SU{2) q which is a g-deformation of the rotation group SU{2). For generic values of q 
the representation theory of SU{2) q is similar to that of SU(2), i.e. every irrep corresponds 
to an integer or half integer spin j = 0, 1/2, . . . and the tensor product of irreps satisfies the 
standard Clebsch-Gordan decomposition. However, if q is a root of unit, q = e 27r4 /( fe+2 ) then 
there is only a finite number of regular irreps corresponding to the spins j = 0, 1/2, . . . , k/2. 
The tensor product of these irreps is a truncated version of the classical CG decomposition, 

(Vn®V^ = &p^> k -^V J (18) 

Vj denotes the vector space of dimension 2j + 1 associated to the irrep with spin j. It 
is interesting to observe that the truncated tensor product (0) satisfies the associativity 
condition (||), namely fT6|| . 

((vi ® v 2 y ® v 3 y = 04 ® (v 2 ® v 3 yy (19) 



This eq. is a consequence of the co-associativity of the comultiplication of the quantum 
group SU(2) q . In more physical terms, eq.(|l^) follows from the non trivial addition rule of 
angular momenta in SU(2) q . The regular irreps have positive g-dimension which is defined 

as 



rj 



(2j+l)/2 _ -(2i+l)/2 

d i = [2J + 1], = g i/2_g-i/2 (2°) 

The q— dimension of an irrep plays a role similar to the eigenvalues w a of the density 
matrix in the sense that irreps with zero g-dimension are thrown away in the tensor product 
just like in the DMRG truncation. Based on this analogy we conjecture that a g-group 
invariant Hamiltonian, like the XXZ open chain with q a root of unit, when studied with 
DMRG methods will yield a density matrix with vanishing eigenvalues corresponding to 
non regular irreps. The DMRG truncation of these states have to agree with the g-group 
truncation of the non regular states ||. 

Quantum groups with q a root of unit are intimately related to rational CFT's (RCFT). 
Indeed in a RCFT there is a finite number of primary fields <f) a { a — 1> • • • , M), which are 
in one-to-one correspondence with the regular irreps of the associated g-group [|16| . Hence 
from the previous relation between the DMRG and g-groups we may expect a relationship 
between RCFT's and the DMRG. More generally in a CFT there are null states in the 
Verma modules of the primary fields, whose norm is zero. As shown by Belavin, Polyakov 
and Zamolodchikov (BPZ) the decoupling of null vectors leads to a set of partial differential 
equations for the conformal blocks of the theory, in terms of which one can construct all 
the correlators of the theory [|17]]. It is tempting to suggest that the BPZ decoupling of 



null vectors is the field theoretical version of the DMRG truncation. On the other hand the 
analogue of the tensor product decomposition is given by the fusion rules of the primary 
fields, 

a <g> <p b = <p c (21) 

where N£ b is an integer which counts how many times the primary field C appears into the 
Operator Product Expansion (OPE) of 4> a and <pb. The associativity of the OPE, i.e. 
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((0a ® 0b) ® 0c) = (0a ® (06 ® 0c)) (22) 

implies a non trivial eq. for the fusion coefficients N£ b namely, 

E N U N L = E N L K ( 23 ) 

a d 

An example of RCFT is given by the SU(2)k WZW model with level k \TB\. The primary 
fields 4>j are labelled by the spin j = 0,1/2, ... , k/2, while the fusion rules are given by the 
eq. fll8D with the translation Vj — > (f)j. Indeed, as shown in reference [16], there is a one-to-one 



correspondence between the SU{2)k WZW model and the quantum group SU(2) q . 

Another aspect of the relation between CFT, Integrable Systems and the DMRG concerns 
the explanation of the exponential decay of the eigenvalues of the density matrix. An 
approach to study this connection is through the relation between the DMRG density matrix 
and the Corner Transfer Matrix (CTM) of Baxter first pointed out by Nishino Jl9|| in his 



application of the DMRG to classical statistical mechanical models in 2D. As shown by 
Baxter |5D| in an integrable system the eigenvalues of the CTM have a very simple structure, 



i.e. they go as a n with n an integer. One can recognize here the exponential decay of the 
eigenvalues of the density matrix [^] , . 

There are still many aspects to clarify in the relation between the DMRG and CFT and 
more generally integrable systems. This could be a fruitful subject in the near future. 



THE DMRG AND THE MATRIX PRODUCT ANSATZS 

The DMRG is a variational method which generates an ansatz for the GS state and 
the excited states. This implies in particular that the DMRG ground state (GS) energy is 
an upper bound of the exact GS energy. The variational ansatz generated by the DMRG 
is of the matrix product (MP) type. This fact was shown by Ostlund and Rommer in 
the thermodynamic limit of the DMRG in the case of the spin 1 chain [7[]. These authors 
proposed that one could get very good results for the GS energy and spin gap by using a 
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MP ansatz which corresponds to a small value of m in the DMRG. The excitations could 
also be constructed as Bloch waves on the MPM state. 

To understand why the DMRG gives rise to a MP state we return to eq. (|T3|) . If ^2) 
and I a) denote basis of the Hilbert spaces associated to B\, B2 and (B1B2)', then the relation 
between these basis is 

mi TU2 

\a) = E u hh,<* \h) 1*2), (a = l,...,m) (24) 

ii=l i 2 =l 

Since B 2 is usually a lattice site • we shall write eq. fl24"|) in the following form, 

\<*)n = E KM M (25) 

where \sn) denotes the local state associated to the site located at the -/V th position of the 
chain, while |a)jv and tv— 1 are the states kept for the blocks of lengths N and N — 1 
respectively. A^Js] is a matrix m x m for each value of s. 

Iterating until reaching the boundary of the chain one gets, 

\a N ,a ) N = (A N [s N ]A N - l [s N - 1 ]...A 1 [s 1 ]) |s x ) . . . | SjV -i} M (26) 

V /a N ,ao 

where the matrix multiplication of the matrices is implicit. |otjv, cxq) is a collection of 

states of an open chain with iV sites labelled by the pair (a/v, «o)- For a closed chain with 
periodic boundary conditions the ansazt becomes 

\iP) N = Tr{A N [s N ) A N ~ l [s N ^} . . . A 1 ^)) \ Sl ) . . . \sn-x) \s N ) (27) 

A further simplication of (p7|) is to assume that all the matrices A n [s n ] are independent 
on n, i.e. A n [s n ] = A[s n ] ( for all n). This assumption can be justified in the thermodynamic 
limit of the DMRG where it reaches a fixed point |7j. However for finite dimensional systems 
and specially for open BC's there will be a non trivial dependence of v4 n [s n ] on n. In this 
sense the DMRG gives a non homogenous MP ansatz. 

In eq.(PB|) we may want that the states \a) form an orthonormal set of states given that 
both \f3) and \s) are orthonormal sets. This implies the following eq. on A^fs] 
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E(<*Mr 45>1 = <W (28) 

As 

which is nothing else than the eq. (|14]) . This eq. expresses the fact that A N relates 
orthonormal basis. But recall that it is not simply a change of basis because we are truncating 
states, i.e. m < m\m2- 

Given the MP ansatzs ( |26| ) and ( p7| ) for open and closed chains respectively we can use 
a standard variational method to find the amplitudes A[s] which minimize the energy of the 



ansatz. In references [0,0 h was shown that when N is large these minimization procedure 
is similar to the one of the DMRG and that in fact there is a hidden density matrix even 
though the algorithm did not try to follow the DMRG method. One way to see this is if one 
define the following transfer matrix, 

^W = E (29) 

s 

T N is a m 2 x m 2 matrix which serves to relate matrix elements of operators between states 
with lengths N and N — 1, for example 

N (a\0\a') N = £ T^, M> n _ x {P\0\P)n-i (30) 

13,13' 



We are assuming in (p0|) that the operator O does not act on the N th site. The normalization 
condition (|28|) implies that T has a right eigenvector with eigenvalue 1 given by 5 a>a ', namely 

12 T aa',f3/3'$/3,l3' = <W' (31) 

It then follows that T has a left eigenvector with eigenvalue 1, i.e. 

X/ Paa'Taa' ,/3/3' = P/3,/3' (32) 



In references [13,14] it was shown that p^ a , can be identified with a density matrix and 
that one is really minimizing the expectation value of the Hamiltonian in the following mixed 

state, 

P N = Y,Paa'\^)N N(a'\ (33) 
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From this point of view the collection of states | a) n of the MP method can be interpreted 
as the most probable ones that contribute to the GS wave function of a system with N + 
1 + N = 2N + 1 sites. 

The important conclusion to be learn from the previous considerations is that the MP 
ansatz leads in a natural way to the DMRG algorithm. This may be interesting regarding 
further generalizations of the DMRG to higher dimensions. 



MATRIX PRODUCTS ANSATZS IN 2D 

The DMRG algorithm can be generalized to ladders (i.e. collections of a finite number 
of chains), and large clusters. This has been done obtaining remarkable results which are 



difficult to obtain with other algorithms |p3]| . However it has also been shown that the 



efficiency of the DMRG disminish with the width of the system [23J. The DMRG algorithm 
appears to be essentially one dimensional in the sense that the RG steps follow a linear 
pattern no matter whether the system is ID or higher dimensional. 

Of course any higher dimensional system can be converted into a ID system by allowing 
non local interactions. However locality seems to be the key of the great performance of the 
DMRG in ID. Hence a truly higher dimensional version of the DMRG should try to keep 
locality as a guideline. That this is in principle possible is suggested by reference |25| where 
a DMRG algorithm is given for a Bethe lattice, whose dimensionality is actually infinite. 



Also recently, Niggemann et al. |26] have constructed a two-dimensional tensor product to 
describe the ground state of a 2D quantum system. Similarly, Nishino and Okunishi have 
proposed a Density Matrix algorithm for 3D classical statistical mechanics models |27| . 



Another approach is suggested by the equivalence of the Matrix Product approach and 
the DMRG for Id or quasi-ld systems. We do not know at the moment what is the formu- 
lation of the DMRG in 2D but we do know that in 2D there are MP states which were first 
constructed by Affleck, Kennedy, Lieb and Tasaki (AKLT) These states are valence 

bond solid states where one connects local spins through local bonds. 
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When trying to generalize the ID MP states to 2D we find that there are two possible 
types of MP states which can be conveniently named as vertex-MP and face-MP states, 
using standard Statistical Mechanics terminology 120 . 



Vertex-Matrix Product ansatzs 

A vertex model in 2D Statistical Mechanics (SM) is a model defined on a square lattice 
and such that the lattice variables i,j, . . . live on the links while the interaction takes place 
on the vertices |2(J . The Bolzmann weight thus depends on 4 variables and the whole 
partition function is obtained by multiplying the Bolztmann weights of all vertices and then 
summing over the values of the lattice variables. The 6 vertex and 8 vertex models are the 
cannonical examples of these types of models which have been shown to be integrable. 

Motivated by these vertex models we shall define a vertex-MP state in terms of a set of 
amplitudes 

Aa*p[s], {a 13, ... = 1, ...,m; s = l,...,m s ) (34) 

where the labels a, f3,j,8 are associated with the links of the square lattice while s labels 
the quantum state, e.g. spin, associated to the vertex where the 4-links a,/3, 7, S meet. 
A„l[s] is a sort of Boltzmann weight of a vertex model. The vertex-MP wave function 
ip(si, S2, ■ ■ ■ , sjv) can be obtained by multiplying all the Boltzmann weights A^'^JsJ and 
contracting and suming over the links variables according to the same pattern of a vertex 
model in Statistical Mechanics [pOj] . Hence the value of the wave function ^(si, S2, ■ ■ ■ , sjv) 
is given by the partition function of a vertex model where the Boltzmann weights depend 
on the value of the local states Sj. This construction for the square lattice is equivalent to 
the so called "vertex-state representation" of Niggemann et al. for the hexagonal lattice [?]. 

This construction resembles the one proposed by Laughlin concerning the Fractional 
Quantum Hall effect (FQHE) pOf . More explicitely, Laughlin proposed in a variational 
wave function i[> m (zi, • • • , z^) for the ground state of the N electrons in the lowest Landau 
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level of a FQHE with filling factor v = 1/m. The norm \i/j m \ 2 of the Laughlin wave function 
can be interpreted as the Boltzman weight of a classical one component plasma constitued 
by N negative charges of magnitude m in a uniform background of positive charges. The 
charge neutrality of the plasma guarantees its stability. 

In our case we also have an associated Statistical Mechanical model given by the Boltz- 
man weights of a vertex model. If we compute the norm of the wave function ip(si, . . . , sjv), 
we can perform the summation over the "spin" indices s, in which case the the norm (ip\ip) 
of the vertex-MP state is given by the partition function of another vertex model whose 
Boltzmann weights are defined as, 

^S' = E^H^?'H (35) 



This R matrix is the 2D version of the T matrix defined in (|29|). The computation of the 
norm of can be in general a difficult task. However, if the model defined by the weights 
(|35|) turns out to be integrable, then we could find the exact norm in the thermodynamic 
limit. 

The face-MP models can be defined in a similar manner by a set of variational parameters 
as in ( P^ ) where now the variables a, . . . are now associated to the vertices of the squares 
while the quantum variable s is associated to the face whose vertices are a,/3,j, 8. This is 
similar to the face or Interaction Round a Face models (IRF) in Statistical Mechanics pO 



Hence in 2D there are two generic ways to produce MP ansatzs which are in fact the 
straightforward generalization of the ID MP ansatzs. These two generalizations suggest to 
use some well know models as the 6- vertex model to test some of the ideas presented above. 

In summary, we have tried to show in this contribution some interesting connections 
among seemingly unrelated methods in condensed matter and field theory. Much remains 
to be done along this direction. 
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